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The  intercalation  of  various  atomic  species,  such  as  hydrogen,  to  the  interface  between  epitaxial 
graphene  (EG)  and  its  SiC  substrate  is  known  to  significantly  influence  the  electronic  properties  of 
the  graphene  overlayers.  Here,  we  use  high-resolution  X-ray  reflectivity  to  investigate  the  structural 
consequences  of  the  hydrogen  intercalation  process  used  in  the  formation  of  quasi-free-standing 
(QFS)  EG/SiC(0001).  We  confirm  that  the  interfacial  layer  is  converted  to  a  layer  structurally  indis¬ 
tinguishable  from  that  of  the  overlying  graphene  layers.  This  newly  formed  graphene  layer 
becomes  decoupled  from  the  SiC  substrate  and,  along  with  the  other  graphene  layers  within  the 
him,  is  vertically  displaced  by  ^2.1  A.  The  number  of  total  carbon  layers  is  conserved  during  the 
process,  and  we  observe  no  other  structural  changes  such  as  interlayer  intercalation  or  expansion  of 
the  graphene  d-spacing.  These  results  clarify  the  under-determined  structure  of  hydrogen  interca¬ 
lated  QFS-EG/SiC(0001)  and  provide  a  precise  model  to  inform  further  fundamental  and  practical 
understanding  of  the  system.  ©  2014  AIP  Publishing  LLC.  [http://dx.doi.org/10.1063/L4899142] 


The  unique  properties  of  graphene  hold  great  potential  for 
use  in  nanoelectronics, ^  and  one  of  the  most  practical  routes 
for  wafer- scale  production  of  high-quality  graphene  is  by  the 
thermal  decomposition  of  SiC  to  produce  epitaxial  graphene 
(EG).5,6  EG/SiC  exhibits  the  electronic  and  structural  proper¬ 
ties  characteristic  of  graphene  and  remains  the  most  promising 
form  of  graphene  for  use  in  state-of-the-art  devices  such  as 
quantum  resistance  standards  and  radio  frequency  field-effect 
transistors.7,8  In  the  case  of  epitaxial  graphene  on  Si- 
terminated  SiC  [EG/SiC(0001)],  however,  the  inferior  electron 
mobility  and  increased  charge  carrier  scattering9-12  are  often 
attributed  to  the  influence  of  the  ubiquitous  (6x/3  x  6x/3)7?30° 
(< 6R3 )  interfacial  layer.9,13-15  Recently,  though,  it  has  been 
demonstrated  that  high-temperature  annealing  of  EG/ 
SiC(0001)  in  H2  leads  to  the  intercalation  of  hydrogen  to  the 
interface.16  In  this  process,  the  hydrogen  breaks  the  Si-C 
bonds  between  the  electronically  inactive,  but  structurally 
graphene-like  interface  layer  and  the  SiC(0001)  (hereafter 
referred  to  as  SiC)  surface,  thereby  decoupling  the  interfacial 
layer  from  the  substrate.  The  resultant  “quasi-free-standing” 
(QFS)  graphene  exhibits  enhanced  carrier  mobility  and  den¬ 
sity,  and  field-effect  transistors  made  from  QFS-EG/SiC  show 
higher  current  saturation  and  transconductance,  and  improved 
extrinsic  current  gain  response.17,18  Furthermore,  these 
improved  properties  are  accompanied  by  a  reversal  in  dopant 
type  from  n-  to  p- type.  This  result  emphasizes  the  considerable 
effect  of  the  local  environment  on  such  low-dimensional  elec¬ 
tronic  stmctures. 

This  resultant  QFS-EG/SiC  structure,  in  which  it  is 
believed  that  the  newly  decoupled  graphene  layer  sits  atop  a 
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H-saturated  hexagonal  SiC  surface  cell,16,18  possesses  an 
interfacial  structure  distinct  from  as-grown  (AG)  EG/SiC. 
The  Si-H  termination  is  thought  to  chemically  and  electroni¬ 
cally  passivate  the  terminal  SiC  layer,  and  the  subsequent 
decoupling  of  the  interfacial  layer  from  the  substrate  elimi¬ 
nates  much  of  its  detrimental  influence.  Regardless,  the  elec¬ 
tronic  properties  of  QFS-EG/SiC  generally  remain  markedly 
inferior  as  compared  to  those  made  from  suspended  gra¬ 
phene,19  Si02-supported  transferred  graphene,20  or  EG 
grown  on  the  carbon-terminated  surface  of  SiC(OOOl).21  It 
therefore  begs  the  question:  What  influence  does  the  newly 
formed  Si-H  terminated  SiC  interface  have  on  the  overlying 
graphene?  Only  with  precise  understanding  of  the  physical, 
chemical,  and  electronic  structure  of  the  interface  can  this 
question  be  addressed. 

In  this  work,  we  investigate  the  atomic  scale  stmcture  of 
H-intercalated  QFS-EG/SiC  using  synchrotron-based  high-re¬ 
solution  X-ray  reflectivity  (XRR).  Our  analysis  reveals 
the  positions  of  the  constituent  atomic  layers  along  the 
SiC[0001]  direction  with  sub- A  resolution  and  verifies  the 
decoupling  and  vertical  displacement  of  the  C-only  interfacial 
layer  from  the  substrate  after  H  intercalation.  Previously, 
some  structural  properties  of  QFS-EG/SiC  have  been  inferred 
using  from  core-level  and  angle-resolved  X-ray  photoelectron 
spectroscopy  (XPS),16  scanning  tunneling  microscopy 
(STM),14  low-energy  electron  microscopy  (LEEM),16  and 
transmission  electron  microscopy  (TEM)  22  However,  these 
approaches  are  either  destructive  or  largely  incapable  of 
accessing  the  buried  interface  structure  with  the  sensitivity 
afforded  by  XRR.  We  also  consider  proposed  alternative 
structural  models,  such  as  the  possibility  of  hydrogen  inter¬ 
calation  between  graphene  layers. 

Two  semi-insulating,  nominally  identical  0.4°  offcut 
6H-SiC  substrates  (II- VI,  Inc.),  cleaved  from  the  same  wafer, 
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were  first  etched  in  80  slm  flowing  H2  at  1520  °C  for  0.5  h  in 
order  to  remove  polishing  damage.  Samples  were  then 
graphitized  under  identical  conditions  (1540  °C  in  lOOmbar 
Ar  for  0.42  h)  in  an  Aixtron/Epigress  VP508  Hot-Wall  CVD 
reactor.  Following  graphitization,  both  samples  were  cooled 
to  1050  °C  in  Ar,  at  which  point  the  pressure  was  increased 
to  900  millibar  of  either  Ar-  or  Pd-purified  H2.  Both  samples 
were  annealed  for  0.5  h  in  these  respective  environments  and 
then  cooled  to  700  °C  before  evacuating  the  chamber. 
Samples  were  removed  from  the  reactor  and  exposed  to  am¬ 
bient  before  cooling  to  room  temperature.  The  in  Ar-  and 
H2-annealed  are  hereafter  referred  to  as  as-grown  AG-EG/ 
SiC  and  QFS-EG/SiC,  respectively.  After  transfer  in  ambi¬ 
ent,  samples  were  then  loaded  into  UHV  and  degassed  for 
lh  at  500  °C  for  low-energy  electron  diffraction  (LEED) 
measurements. 

XRR  measurements  were  performed  under  He  flow  at 
beamline  6-ID-B  at  the  Advanced  Photon  Source  (APS)  at 
Argonne  National  Laboratory.  Measurements  were  performed 
using  =  15.00 keV  (wavelength  X  —  0.8265  A)  X-rays,  follow¬ 
ing  the  methods  described  in  Refs.  23  and  24.  Data  were  fit 
between  qz  =  0A  to  6.6  A-1,  where  qz  =  4nSin(26/2)/X  is  the 
out-of-plane  component  of  the  momentum  transfer  and  26  is 
the  scattering  angle.  Reflectivity  is  plotted  as  a  function  of 
6H-SiC  reciprocal  lattice  units  (SiC  r.l.u)  L  =  (cSic  #z)/27t, 
where  cs ic  =  15.12  A  is  the  6H-SiC  lattice  parameter. 

LEED  patterns  for  control  and  H2-processed  EG/SiC 
[Figures  1(a)  and  1(b)]  verify  the  intercalation  of  the  H  to 
the  interface  as  originally  observed  by  Riedl  et  al}6  Figure 
1(a)  shows  the  typical  pattern  for  as-grown  EG/SiC,  exhibit¬ 
ing  (hk)  spots  indexed  to  the  SiC(10)  and  EG(10)  reciprocal 
lattice  rods,  as  well  as  the  6R3  superstructure  spots  associ¬ 
ated  with  the  interfacial  layer.  After  annealing  in  H2,  the  6R3 
spots  are  suppressed  while  the  EG  peaks  increase  in  relative 
intensity  [Fig.  1(b)] — a  hallmark  that  the  intercalation  pro¬ 
cess  has  occurred.16 

X-ray  reflectivity  data  and  best-fit  models  are  shown  in 
Fig.  2.  The  XRR  data  from  the  AG-EG/SiC  sample  is  compa¬ 
rable  to  that  observed  in  Ref.  13  and  other  similar  works.23,25 
Comparatively,  the  QFS-EG/SiC  shows  orders-of-magnitude 
differences  in  scattered  intensity  at  certain  L  values,  signify¬ 
ing  a  distinct  change  in  the  interfacial  structure.  The  most 
prominent  of  these  differences  are  at  L  ~  5.5  and  L  ~  12.5, 
where  the  QFS-EG/SiC  data  show  deep  minima  where  none 
existed  before  H  intercalation. 

AG-EG/SiC  and  QFS-EG/SiC  electron  density  ( Ne )  pro¬ 
files  derived  from  the  fits  in  Fig.  2(a)  are  shown  in  Figs.  2(b) 
and  2(c),  and  corresponding  interfacial  ball- and- stick  models 
are  presented  in  Figs.  3(a)  and  3(b).  These  profiles  are  calcu¬ 
lated  from  models  comprised  of  k  atomic  layers  defined  by 
their  positions  (zk),  occupancies  (0*),  and  distribution  widths 
(o'*),  as  described  in  Refs.  23  and  26,  and  normalized  to  the 
electron  density  of  an  ideal  Si  layer  within  the  SiC  substrate 
( Ne/NSi  =1).  We  base  this  model  on  the  results  of  our  previ¬ 
ous  work,  in  which  the  interfacial  structure  of  AG-EG/SiC 
was  resolved  using  XRR  in  conjunction  with  X-ray  standing 
wave-enhanced  XPS.13  Using  this  model,  the  data  were  fit 
with  an  interfacial  layer  comprised  of  two  distinct  C  layers 
and  a  five-layer  graphene  film  with  monotonically  reducing 
layer  occupancy  with  increasing  height  above  the  SiC 


FIG.  1.  LEED  patterns  at  130eV  for  AG-  and  QFS-EG/SiC.  SiC(10)  and 
EG(10)  LEED  spots  are  indicated  in  red  and  white,  respectively,  (a)  Before 
intercalation,  the  typical  6R3  superstructure  spots  are  evident,  (b)  After 
intercalation,  the  6R3  spots  are  strongly  suppressed  and  the  EG  spots 
increase  in  relative  intensity. 

surface.  However,  this  initial  model  yielded  very  poor  fits, 
especially  in  the  low-qz  region,  resulting  in  goodness-of-fit 
X2  values27  greater  than  80.  The  source  of  this  original  poor 
fit  is  unclear,  but  the  poor  fits  are  limited  to  the  low-qz 
region,  indicating  the  presence  of  nm-scale,  near-surface 
electron  density  contrast.  We  therefore  employed  a  model 
which  included  a  jellium-like  electron  density  distribution 
within  the  top  few  SiC  bilayers.  This  model  yields  greatly 
improved  XRR  fit  results,  reducing  the  best-fit  j1  value  by 
over  an  order  of  magnitude  to  /2  =  5.27  for  AG-EG/SiC.  The 
electron  density  profile  in  Fig.  2(b)  corresponds  to  the  best- 
fit  AG-EG/SiC  data  in  Fig.  2(a)  and  finds  the  interface  layer 
to  be,  apart  from  the  disordered  component,  statistically 
identical  to  that  found  in  Ref.  13.  Possible  sources  for  this 
disorder  may  be  small  quantities  of  SiOx  inclusions  or 
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L  (SiC  r.I.u) 


FIG.  3.  Ball-and-stick  representations  of  the  conversion  from  AG-EG/SiC  to 
QFS-EG/SiC  as  derived  from  XRR  analysis  in  Figure  2.  H,  C,  and  Si  atoms 
are  colored  in  white,  black,  and  red,  respectively. 

electron  density  profile  is  shown  in  Fig.  2(c).  Best-fit  results 
for  zk,  0£,  and  ok  parameters  for  the  graphene  ( k  =  1-6)  and 
H  layer  ( k  =  H)  are  reported  in  Table  I.  Here,  a  H  layer  whose 
occupancy  was  constrained  to  that  of  the  topmost  Si  layer 
was  added  at  the  nominal  Si-H  bond  distance  of  1.5  A.  Due 
to  its  low  electron  density,  the  fit  results  were  largely  insensi¬ 
tive  to  the  structural  details  of  the  H  monolayer  (ML), 
although  its  inclusion  did  improve  the  best-fit  j1  value  by~l. 
The  total  graphene  coverage  of  the  partially  occupied  layers 
for  QFS-EG/SiC  is  3.0  ML,  which  is  ~1  EG  ML  greater  than 
that  of  the  AG-EG/SiC  sister  sample,  as  would  be  expected 
as  the  interfacial  layer  converts  to  an  undistorted  graphene 
layer.  As  observed  by  comparing  Figs.  2(b)  and  2(c),  the  EG0 
layer  transforms  into  an  EGi  layer  which  is  shifted  upwards 
and  sharpened.  The  characteristic  ~1A  corrugation  of  the 
interfacial  layer13,28  is  reduced  in  amplitude,  making  the  layer 
more  graphene-like,  as  is  consistent  with  STM  studies.14 
Previous  observations  that  Si  and  S2  XPS  components  con¬ 
vert  to  EG  components  corroborate  this  result.16  This  newly 
formed  EGi  is  displaced  from  the  terminal  Si  in  the  SiC  by 
4.2  A,  or  similarly,  from  the  nominal  H  monolayer,  by  2.7  A. 
X2  mapping  of  the  absolute  position  of  the  newly  formed  EGi 
layer  indicates  high  confidence  in  this  value  within  the  limits 
of  the  model:  3-sigma  confidence  levels  for  the  Si-EGx 


FIG.  2.  XRR  data  and  best  fits,  together  with  their  corresponding  electron 
density  profiles,  for  AG-  and  QFS-EG/SiC.  (a)  AG-EG/SiC  (red)  and  QFS- 
EG/SiC  (blue)  data  are  plotted  with  best-fit  curves.  Qualitatively,  the  two 
curves  show  sharp  contrast,  especially  in  the  regions  near  5.5  and 
~12.5,  indicating  a  distinct  change  in  the  interfacial  structure  upon  H  inter¬ 
calation.  (b)  The  electron  density  profile  for  AG-EG/SiC  corresponding  to 
the  fit  in  (a).  A  broad  interfacial  layer  (EG0)  is  comprised  of  two  distinct  C 
distributions  with  an  electron  density  maximum  2.1  A  above  the  terminal  Si. 
(c)  The  QFS-EG/SiC  electron  density  shows  that  the  interfacial  layer  in  (b) 
has  been  converted  to  a  layer  with  graphene-like  density  positioned  4.22  A 
above  the  SiC  surface.  Layer  coverage  and  EG  d-spacing  are  conserved 
upon  H  intercalation.  For  clarity,  the  AG-EG/SiC  data  and  fits  in  (a)  are  off¬ 
set  vertically  by  a  factor  of  103. 

remnant  growth  defects.  Regardless,  this  component  does 
not  appear  to  affect  fit  results  above  the  region  of  L  >  3,  and 
therefore  we  contend  that  the  atomic-scale  structural  analysis 
is  robust.  AG-EG/SiC  fitting  results  for  the  zk,  0^,  and  ok  pa¬ 
rameters  for  the  graphene  {k  =  1-5)  and  interfacial  layers 
( k  =  0)  are  reported  in  Table  I. 

The  XRR  result  for  AG-EG/SiC  was  used  as  the  prelimi¬ 
nary  model  for  the  analysis  of  the  QFS-EG/SiC  data.  The 
best-fit  result  is  shown  in  Fig.  2(a),  and  the  associated 


TABLE  I.  XRR-derived  positions  (zk),  atomic  layer  occupancies  (0/t),  and 
distribution  widths  (o>)  for  AG-  and  QFS-EG/SiC  for  each  k  layer  in  the 
model.  zk  is  reported  with  respect  to  the  topmost  Si  in  the  SiC  substrate  and 
atomic  0^  are  reported  normalized  to  an  ideal  EG/SiC  monolayer;  0  =  1.00 
corresponds  to  38.20  C/nm2  (the  areal  density  of  C  atoms  in  a  graphene 
layer).  Note  that  the  areal  density  of  Si  within  the  SiC  surface  unit  cell  is 
12.17  Si/nm2;  therefore,  a  saturated  Si-H  monolayer  corresponds  to  a  cover¬ 
age  of  0  =  0.32  when  reported  normalized  to  an  ideal  EG/SiC  monolayer. 
Starred  values  were  held  during  fitting. 


k 

AG-EG/SiC 

QFS-EG/SiC 

zk(  A) 

0* 

MA) 

Zk  (A) 

0/t 

<*k  (A) 

H 

1.5* 

0.26* 

0.1* 

0 

f  S2 

1  q, 

2.10 

0.26 

0.1 

2.45 

0.73 

0.40 

1 

5.75 

1.00 

0.15 

4.22 

1.00 

0.1 

2 

9.17 

0.48 

0.11 

7.57 

1.00 

0.14 

3 

12.54 

0.29 

0.1* 

10.95 

0.50 

0.09 

4 

15.84 

0.06 

0.1* 

14.33 

0.25 

0.1* 

5 

19.07 

0.01 

0.1* 

17.64 

0.19 

0.1* 

6 

20.93 

0.08 

0.1* 
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displacement  are  4.22^qo8  A.  In  addition,  the  average  inter¬ 
layer  spacing  between  the  graphene  sheets  for  the  QFS-EG/ 
SiC  EG  is  3.34  A,  which  agrees  closely  with  that  of  AG-EG/ 
SiC  (3.33  A).  This  agreement  verifies  that  the  intercalation 
process  has  essentially  no  effect  on  the  interplanar  spacing  of 
the  EG  film. 

The  2.72+008  H-EGi  displacement  compares  to  within 
0. 1  A  of  recently  calculated  results  using  density  functional 
theory.29,30  Moreover,  recent  work  by  Rajput  et  al.  has 
shown  that  the  Dirac  point  of  monolayer  QFS-EG/SiC  can 
shift  by  nearly  500  meV  as  the  H-EGi  separation  is  varied 
between  2.0  and  3.5  A,  inferring  that  the  hole  density  may 
change  by  orders  of  magnitude  as  a  result  of  subtle  changes 
in  interfacial  displacement.  The  shift  of  the  Dirac  point 
defines  the  Schottky  barrier  height  and  will  determine  the 
practicality  of  employing  the  wide-bandgap  semiconducting 
SiC  substrate  in  a  Schottky  junction  for  application  in,  for 
example,  graphene-based  three-terminal  vertical  field-effect 
transistors.31  The  precise  structural  result  presented  herein  is 
therefore  useful  in  informing  both  computational  and  applied 
efforts  in  determining  the  electronic  influence  of  the  QFS- 
EG/SiC  interface. 

A  final  observation  concerns  the  influence  of  H  intercala¬ 
tion  on  the  interlayer  spacing  of  the  graphene  film  itself.  It 
remains  an  open  question  as  to  whether  the  high  temperature 
processing  of  AG-EG/SiC  in  H2  results  in  the  intercalation  of 
H  species  between  the  graphene  sheets  themselves  in  addition 
to  decoupling  the  interfacial  layer  from  the  substrate. 
Laboratory-based  X-ray  diffraction  and  TEM  results  by 
Tokarczyk  et  al.  suggest  that  the  EG  film  may  expand,  with 
the  observed  p- type  doping  resulting  from  negatively  charged 
H  ions  positioned  between  the  graphene  layers.22  However, 
in  their  diffraction  analysis,  the  authors  do  not  account  for  the 
effect  of  the  surface  truncation  and  interfacial  stmcture  on  the 
X-ray  scattering  profile,  factors  which  we  observe  cause  a 
shift  in  the  absolute  L  position  of  the  graphene  film’s  (0002) 
peak  by  —0.2  SiC  r.l.u  (Fig.  2).  While  our  data  agree  qualita¬ 
tively  with  the  observations  of  Tokarczyk  et  al.,22  by 
accounting  for  the  influence  of  the  interfacial  structure  on  the 
XRR  intensity,  we  find  the  average  interlayer  graphene  spac¬ 
ing  for  AG-  and  QFS-EG/SiC  layers  to  be  statistically  identi¬ 
cal.  We  emphasize  that  while  XRR  is  not  directly  sensitive  to 
the  presence  of  hydrogen  between  the  graphene  sheets,  H 
interlayer  species  are  predicted  to  cause  large  average  gra¬ 
phene  interlayer  expansion  even  for  nanoscale  graphene 
flakes  (EG(0002)  ^3.7-3.9  A)32  or  very  low  H/C  concentra¬ 
tions  (EG(0002)  ^-4.4  A).33  Such  EG  interlayer  expansion  is 
not  observed  in  this  work  despite  the  ~0.1  A  confidence  of 
the  XRR  model  fits  to  the  average  graphene  interlayer  separa¬ 
tion  spacing.  This  suggests  that  if,  indeed,  the  H  intercalation 
process  produces  trapped  interplanar  H  species,  they  do  not 
modify  the  interplanar  graphene  spacing  averaged  over  large 
(~pm)  lateral  length  scales. 

In  summary,  we  used  X-ray  reflectivity  to  compare  the 
structures  of  two  sister  EG/SiC  samples,  one  of  which  under¬ 
went  H  intercalation.  Utilizing  a  previously  established 
structural  model  for  AG-EG/SiC,13  we  are  able  to  verify  the 
physical  decoupling  of  the  interfacial  layer  from  the  SiC  sub¬ 
strate  upon  H-intercalation,  its  displacement  along  the 
SiC [0001]  direction,  and  the  conservation  of  the  total 


coverage  of  overlayer  carbon.  We  report,  with  3-sigma  confi¬ 
dence,  the  vertical  displacement  of  the  newly  formed  EG 
layer  to  be  4.22+ggg  A  above  the  SiC  terminus,  or  2.72+ggg 
A  above  the  nominal  H  layer.  Our  structure  is  in  close  agree¬ 
ment  with  recent  calculations29  and  has  implications  on  the 
electronic  behavior  of  the  graphene  overlayers  themselves. 
This  investigation  of  the  structural  consequences  of  H  inter¬ 
calation  of  EG/SiC  precisely  determines  key  physical  param¬ 
eters  (e.g.,  SiC-EG  displacement  and  interlayer  graphene 
spacing)  that  are  thought  to  critically  influence  technologi¬ 
cally  relevant  properties  such  as  Dirac  point  shift  and 
Schottky  barrier  height.  Furthermore,  this  methodology 
should  prove  valuable  in  the  continued  exploration  of  the 
myriad  of  atomic  species  that  can  be  intercalated  into  EG/ 
SiC  interface  (see  Ref.  34  and  references  therein),  many  of 
which  exhibit  unique  influence  on  the  electronic  properties 
of  the  overlying  epitaxial  graphene.  Indeed,  tailored  interfa¬ 
cial  design  via  intercalation  may  ultimately  prove  to  be  an 
important  tool  in  the  engineering  of  properties  of  epitaxial 
graphene  on  SiC.  Similar  studies  on  such  analogous  systems 
will  prove  critical  to  the  development,  understanding,  and 
utilization  in  these  promising  new  graphene  systems. 
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